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ABSTRACT 

An extensive program for the calculation of galactic cosmic-ray propagation has been developed. 
Primary and secondary nucleons, primary and secondary electrons, and secondary positrons are in¬ 
cluded. The basic spatial propagation mechanisms are (momentum-dependent) diffusion, convection, 
while in momentum space energy loss and diffusive reacceleration are treated. Fragmentation and 
energy losses are computed using realistic distributions for the interstellar gas and radiation fields. 

INTRODUCTION 

The main motivation for developing this code (Strong and Moskalenko 1997) is the prediction of dif¬ 
fuse Galactic gamma rays for comparison with data from the CGRO instruments EGRET, COMPTEL 
and OSSE. This is a development of the work described in Strong and Youssefi (1995). More generally 
the idea is to develop a model which self-consistently reproduces observational data of many kinds 
related to cosmic-ray origin and propagation: direct measurements of nuclei, electrons and positrons, 
gamma rays, and synchrotron radiation. These data provide many independent constraints on any 
model and our approach is able to take advantage of this since it must be consistent with all types of 
observation. We emphasize also the use of realistic astrophysical input (e.g. for the gas distribution) 
as well as theoretical developments (e.g. reacceleration). The code is sufficiently general that new 
physical effects can be introduced as required. The basic procedure is first to obtain a set of propaga¬ 
tion parameters which reproduce the cosmic ray B/C ratio, and the spectrum of secondary positrons; 
the same propagation conditions are then applied to primary electrons. Gamma-ray and synchrotron 
emission are then evaluated. Models both with and without reacceleration are considered. 

DESCRIPTION OF THE MODEL 

The models are three dimensional with cylindrical symmetry in the Galaxy, the basic coordinates 
being ( R,z,p ) where R is Galactocentric radius, z is the distance from the Galactic plane and p is 
the total particle momentum. The numerical solution of the transport equation is based on a Crank- 
Nicholson implicit second-order scheme. In the models the propagation region is bounded by z — Zh 
beyond which free escape is assumed. A value Zh = 3 kpc has been adopted since this is within the 
range which is consistent with studies of 10 Be/Be and synchrotron radiation. For a given zh the dif¬ 
fusion coefficient as a function of momentum is determined by B/C for the case of no reacceleration; 
with reacceleration on the other hand it is the reacceleration strength (related to the Alfven speed va) 
which is determined by B/C. Reacceleration provides a natural mechanism to reproduce the B/C ratio 
without an ad-hoc form for the diffusion coefficient (e.g., Heinbach and Simon 1995, Seo and Ptuskin 
1994). The spatial diffusion coefficient for the case without reacceleration is D = p/Jo below rigidity 
Po, PA)(p/po ) 8 above rigidity po, where p is the particle speed. The spatial diffusion coefficent with 
reacceleration assumes a Kolmogorov spectrum of weak MHD turbulence so D = (iC>o(p/po ) 8 with 
8=1/3 for all rigidities. For this case the momentum-space diffusion coefficient is related to the spa¬ 
tial coefficient (Seo and Ptuskin 1994). The injection spectrum of nucleons is assumed to be a power 
law in momentum for the different species, dq(p)/dp p~^ for the injected density, corresponding 
to an injected flux dF(E k )/dE k p Y where E k is the kinetic energy per nucleon. 

The entire calculation is performed with momentum as the kinematic variable. The interstellar 


fraction of the gas is taken as 0.11 by number. The interstellar radiation field for inverse Compton 
losses is based on stellar population models and IRAS and COBE data, plus the cosmic microwave 
background. The magnetic field is assumed to have the form 6e^ z l/ 5kpc )~^/ 20kpc ) pG. Energy losses 
for electrons by ionization, Coulomb, bremsstrahlung, inverse Compton and synchrotron are included, 
and for nucleons by ionization and Coulomb interactions. The distribution of cosmic-ray sources is 
chosen to reproduce the cosmic-ray distribution determined by analysis of EGRET gamma-ray data 
(Strong and Mattox 1996). The bremsstrahlung and inverse Compton gamma rays are computed self- 
consistently from the gas and radiation fields used for the propagation. The 7t°-decay gamma rays 
are calculated directly from the proton and Helium spectra using the method of Dermer (1986). The 
secondary nucleon and secondary e ± source functions are computed from the propagated primary 
distribution and the gas distribution, and the anisotropic distribution of in the system was taken 
into account. 

ILLUSTRATIVE RESULTS 
Some results obtained are shown in the Figures. 

Fig. 1: Secondary nucleons. The energy depen¬ 
dence of the B/C ratio can be reproduced with 
Dq = 2.0xl0 28 cm 2 s -1 , 8 = 0.6, Po = 3 GV/c with¬ 
out reacceleration (thick line) and Do = 4.2 x 10 28 
cm 2 s -1 , va = 20 km s 1 with reacceleration (thin 
line). Dashed lines are modulated to 500 MV. Data 
from Webber et al. (1996). 
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Fig. 2: Left panel: the local proton spectrum for injection index 2.15 (thin solid line), 2.25 (thick 
solid line) without and with reacceleration respectively, compared with the measured ‘interstellar’ 
spectrum (dashed line: Seo et al. 1991, dashed-dot line: Mori 1997). Right panel: the Helium spec¬ 
trum with injection index 2.25 (thin solid line), 2.45 (thick solid line) without and with reacceleration 
respectively, compared with the measured ‘interstellar’ spectrum (dashed line: Seo et al. 1991, 
dashed-dot line: Engelmann et al. 1990). The spectra are well reproduced up to 100 GeV. 


Secondary positrons and electrons: Using the primary proton and Helium spectra, the propaga¬ 
tion of secondary e^ has been computed. The e + spectrum (Fig 3) and the positron fraction (Fig 4) 
agrees well with the most recent data compilation for 0.1-10 GeV; for more details see Moskalenko 
and Strong (1997). 

Primary electrons: The spectrum of primary electrons is also shown in Fig 3. The adopted elec¬ 
tron injection spectrum has a power law index -2.1 up to 10 GeV; this is chosen using the constraints 





around 10 GeV where solar modulation is small and it also satisfies the constraints from J. 

e + +e 

Above 10 GeV a break is required to at least -2.4 for agreement with direct measurements (which 
may however not be necessary if local sources dominate the directly measured high-energy electron 
spectrum). 
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Fig. 3: Spectra of secondary positrons and electrons, and of primary electrons. Full lines: our 
model with no reacceleration (left panel) and with reacceleration (right panel). Dash-dotted lines 
(Protheroe 1982): lower is his leaky-box prediction for e + , upper is his electron spectrum. 
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Fig. 4: Positron fraction for model with no reacceleration (left panel) and with reacceleration (right 
panel). Dashed and dash-dotted lines: Protheroe (1982) predictions of the leaky-box and diffusive 
halo models respectively. The data collection is taken from Barwick et al. (1997). 


The synchrotron spectrum (Fig 5) at high Galactic latitudes is important since its shape con¬ 
strains the shape of the 1-10 GeV electron spectrum. An injection index -2.1 (without reacceleration) 
is the steepest which is allowed by the radio data over the range 38 to 1420 MHz. As illustrated, an 
index -2.4 as often used (e.g. Strong et al. 1996) gives a synchrotron spectrum which is too steep. 

The modelled ganuna-ray spectrum for the inner Galaxy, illustrated here for the case of no 
reacceleration (Fig 6), fits well the COMPTEL and EGRET data up to 500 MeV beyond which there 
is the well-known excess not accounted for by 7t°-dccay with the standard nucleon spectrum. The 
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Fig. 5: The synchrotron spectrum at the NGP, 
compared to predictions for electron injection in¬ 
dices -2.1 (dashed line) and-2.4 (solid line). 


Fig. 6: The y-ray spectrum for the inner Galaxy, 
330° < / < 30°, \b\ < 5°. EGRET data: Strong 
and Mattox (1996). 




Fig. 7: The model y-ray longitude profile (\h\ < 5°) for bremsstrahlung (histogram), inverse Compton 
(middle cwve) and total (upper curve), for energies 1-3, 3-10 and 10-30 MeV. 


COMPTEL spectrum from a recent analysis by Strong et al. (1997) is both lower and flatter than 
previously published. It fits well to the flatter electron injection spectrum required by synchrotron 
data (see above). Inverse Compton is dominant below 10 MeV, bremsstrahlung becomes important 
for 3-200 MeV. The lower bremsstrahlung combined with 7 t°-dccay leads to a good fit to the steep rise 
in this range, in contrast to previous attempts to model the spectrum with a steeper bremsstrahlung 
spectrum (Strong et al. 1996). Below 10 MeV inverse Compton becomes increasingly dominant, 
and bremsstrahlung is negligible below 1 MeV. The longitude profile (Fig 7) at low latitudes from 
this model is compared with that from COMPTEL elsewhere (Strong et al. 1997). The large energy 
dependence of the relative importance of inverse Compton and bremsstrahlung is evident. 

More details can be found on http://www.gamma.mpe-garching.rnpg.de/~aws/aws.html 
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